##################################################
###Create expression data for hugues
##################################################

load("/data4/bsi/neuro/s102384.spore4/DASL_Jenkins/eqtl/eqtl_all.RData")
snps <- read.table("/data4/bsi/neuro/s102384.spore4/DASL_Jenkins/eqtl/eqtl.raw",sep=" ",header=T)
ids <- snps$FID

hugues <- eqtl.t[eqtl.t$ID %in% ids,]
write.table(hugues,"/data4/bsi/neuro/s102384.spore4/DASL_Jenkins/eqtl/hugues/eqtl_expression.txt",sep=" ",row.names=F,col.names=T,quote=F)
